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. We examine the properties of relic neutralinos in the regions of the parameter 

space selected by recent fits to all the electroweak observables within the Min- 
imal Supersymmetric Standard Model. We discuss the relic neutralino cosmic 
abundance and the direct detection rates for these most likely supersymmet- 
ric configurations. We employ the relevant experimental bounds to constrain 
r S^ . the set of the best-fit configurations and discuss the discovery potential of the 

direct search for relic neutralinos. 
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I. INTRODUCTION 



The Minimal Supersymmetric extension of the Standard Model (MSSM) has been anal- 
ysed by many authors in a variety of different theoretical schemes. These various scenarios 
range from the most model-independent one, whose parameters at the Mz scale are con- 
strained only by experimental bounds, with no assumptions about properties at higher scales, 
to more sophisticated ones, where the quantities at Mz are derived from a few parameters 
at the Grand Unification (GU) scale (Mqut — O(10 16 ) GeV), with the further requirement 
that the ElectroWeak Symmetry Breaking (EWSB) occurs radiatively. 

The specific hypotheses that are assumed at Mqut are very crucial for the phenomeno- 
logical properties of the supersymmetric particles at the Mz scale. One often assumes that 
the soft-breaking mass parameters (scalar masses, gaugino masses and trilinear couplings) 
unify at Mqut as well as the gauge couplings and the b-r Yukawa couplings. It is worth 
remarking that relaxing some of these unification requirements may modify substantially the 
supersymmetric phenomenology . The implications of a relaxation of the soft scalar mass 
unification at Mqut on the properties of relic neutralinos (relic abundance and detection 
rates) are discussed in Refs. 

Unfortunately, at present we do not have enough information from experimental data on 
which would be the most reliable supersymmetric scheme. In view of this situation, partic- 
ularly attractive is the approach developed recently, where the parameters of the MSSM at 
Mz are considered as independent variables in a global fit to the electroweak data without 
any further theoretical assumption |5HlO[. The fit to all the available experimental data 
singles out some specific regions of the supersymmetric parameter space as the most likely 
ones. Furthermore, this analysis, which is performed at the Mz scale, may provide useful 
insight into the most plausible theoretical scenario at the M GUT scale. 

In the present note we examine which phenomenological consequences this global fit 
entails, within the MSSM, for relic neutralinos. For the most likely domains of the parameter 
space we discuss the features of the relevant neutralino relic abundance and of the event 
rates for the direct search. 



II. x 2 ANALYSIS WITHIN THE MSSM 

In the global fits of the electroweak data performed in the framework of the MSSM a key 
role is played by the ratio Rb = T(Z — > bb) /T(Z — > hadrons) whose experimental value shows 
a discrepancy (at the level of 3.2 standard deviations) with respect to the value predicted by 
the Standard Model (SM). In fact, the experimental result R e b xp = 0.2211 ± 0.0016 has 
to be confronted to the SM prediction Rijj M = 0.2160 for m t = 170 GeV (or R$ M = 0.2158 
for m t = 180 GeV). The MSSM may resolve this disagreement (at least partially), since the 
contributions of some supersymmetric loop diagrams add positively to the SM prediction 
T2],|6|-[T0| . More specifically, substantial supersymmetric contributions occur in two regions 



of the supersymmetric parameter space: (i) for small tan/5, when both the chargino and one 
top squark are light, and (ii) for large tan/3, when the CP-odd Higgs boson A is light and/or 
both the chargino and one top squark are light (also light neutralino and bottom squark 
would help in increasing Rb). As usual tan/3 is defined as tan/3 = t^/fi, where v\ and v<i are 
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the vacuum expectation values of the Higgs doublets Hi and if 2, which give masses to the 
down-type and up-type quarks, respectively. 

Thus, including Rb in the overall set of the E W data makes the x 2 fit within the MSSM 
superior to the one within the SM. At the same time, the experimental range for Rb, com- 
bined with other important inputs (b — > 57, bound on rji/,, etc.) is very instrumental in 
placing stringent bounds on the supersymmetric parameter space. Due to the properties 
recalled above, it is natural that the best fit of the electroweak data within the MSSM 
automatically selects either (i) very small values of tan/3 (i.e. close to the quasi-IR fixed 
point for a given value of m t ) or (ii) very large values of tan/3 (i.e. of order m t jm^) with 
appropriate ranges for the masses of A boson, chargino, top squark, neutralino and bottom 
squark 0-gO|. 

It is worth recalling that the fit within the MSSM also has the advantage of providing 
a value for the strong coupling constant ct s (Mz), which is lower than the one found with a 
SM fit, and in agreement with the measurements of a s at q 2 <C M| [P^-TT3|]. Actually, this 
is expected, since lower values of a s (Mz) and additional contributions to the Z hadronic 
width are related properties J/J-^[L6|. 

A few words of caution are in order here. First, also the experimental value of the ratio 
R c = T(Z -> cc)/T(Z -> hadrons): R e c xp = 0.1598 ± 0.0070 |TJ shows a disagreement 
(at the level of 1.7 standard deviations) with the SM prediction: R^ = 0.172. Most 
disturbingly, this disagreement in R c cannot be settled within the MSSM. These problems 
affecting Rb and R c could be due to systematic errors in the experimental measurements, 
and might therefore disappear after the analysis of new data. However, it has to be pointed 
out that, even if R c is constrained to the SM value, still the measured value of Rb shows a 



disagreement of 2.6 standard deviations with the SM value [11]. In the present note we focus 



our attention on the implications that the Rb problem would have for relic neutralinos, if the 
measured effect were real and could be accounted for, at least partially, by supersymmetric 
effects within the MSSM. 

As was mentioned above, an important ingredient for a good x 2 & t small tan/? is provided 
by a positive supersymmetric contribution to Rb, when the chargino is light. This feature 
is somewhat weakened |T7| by the new LEP 1.5 lower bound on the chargino mass m x ± > 
65 GeV (if m x ± — m x > 10 GeV, m x denotes the neutralino mass) fllSj . However, this new 



limit does not modify appreciably the conclusions of the x 2 analysis of Ref. pi, which we 
are going to use in the present note, since the lowest values of x 2 at small tan/3 occur at 
larger values of the chargino mass, m x ± ~ (75-90) GeV [ID . 



As far as large values of tan/3 are concerned, one should recall that significant constraints 
may potentially follow from the Z — > bbA process [113. This argument was used in Ref. 



to conclude that large values of tan/3 are disallowed by present experimental data. However, 
much larger statistics and/or more refined analyses of the LEP experimental data on four 
6-jets are still required, before solid conclusions on this point may be drawn (see also Ref. 
[T0|1 ) . Here we consider the large-tan/? sector of the parameter space as a region deserving 
much attention, because of the indications coming from the fit within the MSSM and in view 
of the very interesting perspectives that it offers for further investigation both at accelerators 
and with other experimental means (direct and indirect searches for relic neutralinos). 

In the analysis presented in this note we take as representative values for case (i), i.e. 
small tan/3, and case (ii), i.e. large tan/?, the values: tan/3 = 1.4 and tan/3 = 50, respectively. 
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Now we briefly mention the main features of the best fit of Ref. ||, which are relevant 



to our subsequent discussion on the neutralino phenomenology [^TJ. Here all Z partial- 



widths and asymmetries at the Z pole are calculated and fitted using as free parameters 
the Higgs- mixing parameter fi, the ratio M 2 //i (M 2 is the SU(2) gaugino mass), the mass of 
the CP-odd Higgs boson tua, the three soft scalar masses for the third family, the soft top 
trilinear coupling A t and the strong gauge coupling a s . The usual relation between U(l) 
and SU(2) gaugino masses M\ = | tan 2 6 , vt/M 2 is assumed. For details concerning the best-fit 
procedure we refer to [j^JTOfl . Various values for the top quark mass were considered in the 



fit of Ref. ||. In the present note we set it to m t = 170 GeV. 

For small tan/? (tan/5 = 1.4) the best fit (x 2 — 12-13) does not occur for the smallest 
possible value of the chargino mass, but rather for the range m x ± ~ (75-90) GeV. In fact, for 
such chargino masses the ratio Rb obtains its largest value, Rb ~ 0.2180, from the properties 
of the couplings occurring in the supersymmetric loop corrections to Rb, for one sign of fi, as 
explained in Ref. [10 1. In the configurations providing the best fit (x 2 — 12-13), the lighter 



stop is very light indeed: Mjr 2 ~ 50 GeV, the neutralino mass turns out to be constrained 
in the range m x ~ (25-40) GeV, and the mass of the Higgs boson A is larger than about 
200 GeV. If we consider all the configurations that satisfy the condition x 2 — 16, then 
the chargino and neutralino masses are in the ranges: 70 GeV < m x ± < 200 GeV and 25 
GeV < m x < 100 GeV, respectively. The neutralino may have quite different compositions, 
from very pure gaugino-like to the other extreme of very pure higgsino-like. In the next 
section we will discuss the neutralino relic density. However, we anticipate here which are 
the modifications to the features of the best-fit previously discussed, when this best-fit 
analysis is implemented with the following constraints: (a) the neutralino is the Lightest 
Supersymmetric Particle (LSP), (b) the cosmological bound Qh 2 < 1 is satisfied. Condition 
(a) is imposed here by requiring the sfermion masses to be larger than the neutralino mass 
m x by at least 10%, i.e. mj > 1.1 m x . This requirement prevents some propagators in 
the elastic cross sections from becoming singular and also avoids possible complications 



due to co-annihilation [22|. It turns out that condition (b) is very effective, particularly 
in disallowing many configurations, since, at small tan/3, Q x h 2 is usually very large. Thus, 
adding the two conditions (a) and (b) somewhat worsens the best-fit quality: the minimal x 2 
becomes 13.5 at m x ± ~ 90 GeV and rriA — 100 GeV. In that region of the parameter space 
m x ~ 40 GeV. This last value is expected, because the cosmological bound is satisfied due 
to the very pronounced enhancement of the neutralino pair annihilation cross-section at the 
Z-pole (m x ~ rriz/2). The region of the parameter space satisfying the condition x 2 < 16 
has chargino and neutralino masses constrained in the ranges 65 GeV < m x ± < 110 GeV 
and 25 GeV < m y < 80 GeV. 

Let us now turn to the large-tan/? case (tan/? = 50). The best fit (x 2 — 12-13) occurs for 
the smallest possible values of the chargino mass and of the A-boson mass: m x ± ~ 65 GeV, 
rriA — 55 GeV and for a light neutralino, m x ~ 50 GeV. This allows large values for R b , 
up to Rb = 0.2181. Let us stress that the relevant chargino and neutralino compositions 
are mainly higgsino-like, as implied by the mass relation m x ~ m x ±. If we consider all the 
configurations that satisfy the condition x 2 < 16, then the masses of chargino and neutralino 
are in the ranges: 65 GeV < m v ± < 250 GeV, 50 GeV < m v < 190 GeV. The cosmological 
condition Qh 2 < 1 is not effective in constraining the parameter space here, as is expected 
at large values of tan/3. 
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III. RELIC NEUTRALINOS 



We now discuss the properties of the relic neutralinos, should the supersymmetric pa- 
rameter space be the one singled out by the previously described x 2 fit- I* 1 particular we 
discuss the neutralino relic abundance and detection rates for the direct search, where a relic 
neutralino may be detected by the energy released by its elastic scattering off a nucleus of 
an appropriate set-up [0-|25|]. 

Some words of caution are in order here, concerning the evaluation of the direct detection 
rates. It is well known that all the calculations for detection rates, in addition to the uncer- 
tainties due to the supersymmetric scheme, are also affected by many uncertainties in some 
astrophysical and cosmological properties. This is the case for the neutralino velocity distri- 
bution in the halo (assumed here to be Maxwellian with a r.m.s. velocity of 270 ±24 km ■ s _1 



I, for its escape velocity (v esc = 650 ± 200 km • s p7|) and for the velocity of the Sun 



around the galactic centre (v Q = 232 ± 20 km • s _1 |2B|). All these velocities are affected 
by large error bars, which in turn introduce significant uncertainties in the detection rates, 
especially for light neutralinos. Furthermore, an even more serious uncertainty concerns the 
value of the local (solar neighbourhood) dark matter density p\. A recent determination 
of pi, based on a flattened dark matter distribution and microlensing data, gives the range 
Pi = 0.5lloj7 GeV • cm~ 3 P5[ . This introduces a central value significantly larger than 
previous determinations, for instance the one of Ref. [29]: pi = 0.3 ±0.1 GeV • cm -3 , and 



then entails larger detection rates. Notice that the error bars affecting pi are quite large. 

Once a specific value for the local density of the total dark matter p\ is picked up, we are 
still confronted with the problem of assigning a value to the neutralino local density p x . Here, 
to determine the value of p x to be used in the detection rates, we adopt the following rescaling 
recipe [|30|: for each point of the parameter space, we take into account the relevant value 
of the cosmological neutralino relic density. When Sl x h 2 is larger than a minimal (Qh 2 ) min , 
compatible with observational data and with large-scale structure calculations, we simply 
put p x = pi. When fl x h 2 turns out to be less than (Qh 2 ) m i n , and then the neutralino may 
only provide a fractional contribution Q x h 2 / (Qh 2 ) min = £ to Qh 2 , we take p x = p£. The 
value to be assigned to (Qh 2 ) min is somewhat arbitrary, in the range 0.03 < (Qh 2 ) min < 0.2. 

In the present paper, as far as the values of the aforementioned astrophysical and cos- 
mological parameters are concerned, we use the two sets of values reported in Table 1. Set I 
corresponds essentially to the central values for the various parameters, whereas set II corre- 
sponds to those values of the parameters, which, within the relevant allowed ranges, provide 
the lowest detection rates (once the supersymmetric variables are fixed). Set I is employed 
to discuss the average predicted rates, whereas set II is used to derive constraints on super- 
symmetric configurations, whenever the theoretical rates exceed the present experimental 
bounds. 

As a further source of uncertainties affecting the detection rates, we finally mention the 
strength of the Higgs-nucleon coupling, which plays an important role in the neutralino- 
nucleus scattering. Here, for the relevant coupling we adopt the determination reported in 
Ref. g. 

Many experiments for direct detection of massive weak-interacting dark matter particles 
are under way, using various materials f25" |. In the present note we consider detection rates 
for germanium detectors |35f| . The most consistent procedure for comparing predicted 
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TABLE I. Values of the astrophysical and cosmological parameters relevant to direct detection 
rates. V r .m. s . denotes the root mean square velocity of the neutralino Maxwellian velocity distribu- 
tion in the halo; V esc is the neutralino escape velocity and Vq is the velocity of the Sun around the 
galactic centre; pi oc denotes the local dark matter density and (Qh 2 ) m i n the minimal value of Qh 2 . 
The values of set I are the median values of the various parameters, the values of set II are the 
extreme values of the parameters which, within the physical ranges, provide the lowest estimates 
of the direct rates (once the supersymmetric parameters are fixed). 





Set I 


Set II 


V r . m . s (km ■ s" 1 ) 


270 


245 


F esc (km • s _1 ) 


650 


450 


^(km-s- 1 ) 


232 


212 


p; oc (GeV • cm" 3 ) 


0.5 


0.2 


{yih ) m in 


0.1 


0.2 



detection rates with experimental data is to consider the differential rate dRdirect/ dE ee over 
the whole range of the electron-equivalent energy E ee (or, alternatively, of the nuclear re- 
coil energy). Here, to simplify the discussion, we present our results only in terms of an 
appropriate integrated rate Rdirect- On the basis of the experimental spectra of Ref. |35j 
we derive, as the most stringent upper bound (at 90% C.L.) for light neutralinos (20 GeV 



< m x < 50 GeV), the value p| 



c8 keV 



ra Key 

\ dE ee dR e d ZJdE ee < 0.9 events/(kg • day). (1) 

J6 keV 

In principle, for heavier neutralinos (m x larger than 50 GeV), which are expected to produce 
much flatter energy spectra, a more stringent bound (at 90% C.L., still using the data of 
Ref. [p5p is provided by 

rU keV 

/ dE ee dR e d Zct/dEee < 0.1 events/(kg ■ day). (2) 

J13 keV 

However, since the limit of Eq. @ introduces only a marginal improvement over that of 
Eq. ([]]), for simplicity, in the present paper we only use the bound of Eq. ([]]), which we 
denote by R d ^ ect , and compare it with the relevant theoretical integrated rate Rdirect = 

I& keV dE ee dRdirect/ dE ee . 



IV. RESULTS AND CONCLUSIONS 

Let us now turn to the presentation of our results, based on the supersymmetric con- 
figurations selected by the global fit to the electroweak data of Refs. |8|jll|. We start from 



the small-tan/5 case. In Fig. 1 we show the neutralino relic abundance Q x h 2 versus Rdirect, 
for tan/5 = 1.4. As was mentioned in Sect. 2, most of the configurations singled out by 
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the best-fit procedure at this value of tan/? provide large values of Q x h 2 , in excess of the 
cosmological bound. The scatter plot shown in Fig. 1 displays only those supersymmetric 
configurations that satisfy the condition Qh 2 < 1. The rate Rdirect reported in this figure 
was evaluated, using set I of Table 1. The vertical line corresponds to the experimental 
bound defined in Eq. (|1|). 

We notice in Fig. 1 that many neutralino configurations provide a value for the relic 
abundance within the optimal range for Cold Dark Matter (CDM), 0.1 < QcDMh 2 < 0.3, 
which is suggested by a number of observational data f3~Tfl . Many other supersymmetric 
compositions provide an even larger contribution to the neutralino relic abundance. As 
expected on general grounds, the predicted direct rates are all below the present experimental 
bound of Eq. ([]]). However, it is interesting to note that, by improving the experimental 
sensitivity of the direct experiments by about one order of magnitude (which appears to 
be within reach in the near future), one can start the investigation of a sizeable number of 
neutralino configurations (including many of cosmological interest). 

Let us now turn to the large-tan/3 case. In Fig. 2 we display our results for the neutralino 
relic abundance versus the neutralino mass, as a scatter plot over all the configurations of 
the best-fit analysis with \ 2 — 16, for tan/3 = 50. In this scatter plot the points turn out 
to be grouped in some characteristic branches. This feature is due to the grid used in the 
X 2 fit: along each branch the ratio M 2 /\/j,\ has a fixed value, and each branch stops when 
m x ± reaches the maximal value considered in the fit: m x ± = 200 GeV {m x ± = 250 GeV for 
the value M 2 /\n\ = 1). The highest branch corresponds to the value M 2 /|/x| = 0.5 (gaug- 
ino dominance in the neutralino composition), the intermediate one belongs to the value 
M 2 / 1 /i I = 1 (large gaugino-higgsino mixing), and the lowest branch contains configurations 
with M 2 j\\i\ = 3, 10 (higgsino dominance). The different nature of the neutralino composi- 
tions along the different branches explains the significant differences in their £l x h 2 values. As 
expected, only gaugino-like neutralinos may have values of relic abundance in the desirable 
range 0.1 < QcDAih 2 < 0.3. Actually, only very few configurations within the total best-fit 
sample fall into this range. However, it is rewarding that a large number of states provide 
a relic abundance of some cosmological interest, i.e. in the range Q x h 2 > 0.025. 

Let us now discuss the event rates for direct detection of relic neutralinos. Figure 3 
displays two scatter plots of the rate Rdirect as a function of m x . The rate Rdirect shown in 
section (a) was evaluated using set I in Table 1, the one shown in section (b) was evaluated 
using set II. In Fig. 3 configurations of different cosmological relevance are denoted by 
different symbols. The shape of the rate Rdirect in these plots is reminiscent of the shape 
of the relic abundance displayed in Fig. 2, the reason being that the rescaling of the local 
neutralino density is very effective here. Of the two sets of points which branch off at m x ~ 80 
GeV, the upper branch is mainly a combination of the configurations with M 2 j\[i\ = 0.5, 1, 3, 
whereas the lower branch contains configurations with M 2 j\[i\ = 10. In fact, the neutralino- 
nucleus cross section that enters in the detection rate has the effect of somewhat grouping 
together the configurations belonging to the three smallest values of M 2 /\fi\ and of further 
suppressing the rates for the purest higgsino states with M 2 /\fi\ = 10. Of the two plots 
shown in Fig. 3, the one of section (a), which was obtained by using the values of set I for 
the astrophysical parameters, is indicative of what would be expected from direct detection, 
but suffers from uncertainties due to the large error bars affecting the various parameters. 
The plot of section (b) is the one which is relevant, when we address the question whether the 
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present experimental bound R^ect already allows us to set some constraints on our sample 
of best-fit supersymmetric configurations. As a matter of fact, the scatter plot displayed in 
Fig. 3b represents a conservative estimate of the predicted rates. Thus we may conclude 
that the neutralino configurations, whose rates are above R e ^ ect in the plot of section (b), 
are excluded by the present experimental bound from direct detection experiments |[38| . The 
number of the disallowed configurations is almost one tenth of the total number of about 
3,000 configurations. 

In Fig. 4 we display, in an R^ versus x 2 plot, the frontier of the region that contains the 
whole sample of the best-fit configurations of Ref. (this frontier is denoted by a solid line). 
The dashed curve delimits the region where are located those configurations of the original 
sample that are excluded by R^rect on the basis of the results shown in Fig. 3b. However, 
we stress that the region inside the dashed curve is still populated by a large number of 
other allowed configurations. 

Some noticeable properties emerge from our results: 

• The configurations excluded by the present upper bound Refect are n °t located in the 
minimum of x 2 ( see Fig- 4). Actually, from details of our analysis which are not dis- 
played in Fig. 4, it turns out that most of the excluded configurations are concentrated 
in the range 15 < x 2 ^ 16. The configurations that give the minimal values of x 2 have 
a very pure higgsino composition, and then provide a small detection rate. Therefore, 
the constraint introduced by R^rect^ although instrumental in limiting the parameter 
space, does not nevertheless modify the good quality of the fit significantly. 

• Even a modest improvement in the sensitivity of the direct search experiments allows 
the exploration of a large number of neutralino configurations. This is clear from 
Fig. 3b, which shows that the Rdirect vs. rn x plot is highly populated around the line 
denoting the present upper bound R e J^ ect - For example, our calculations predict that, 
once the experimental sensitivity is improved by an order of magnitude with respect to 
the present one, more than 60% of the total number of the best-fit configurations can 
be explored. Some of them belong to the region of minimal x 2 ■ Thus, the discovery 
potential of the direct search for relic neutralinos is quite remarkable. 



Acknowledgements 

We thank Stefan Pokorski, Piotr Chankowski and Nicolao Fornengo for interesting 
discussions. This work was supported in part by the Research Funds of the Ministero 
dell'Universita e della Ricerca Scientifica e Tecnologica. A fellowship of the Universita di 
Torino is gratefully acknowledged by G. Mignola. M. Olechowski would like to acknowledge 
the support of the Polish Committee of Scientific Research. 



8 



REFERENCES 



[1] M. Olechowski and S. Pokorski, Phys. Lett. B344 (1965) 201. 

[2] N. Polonsky and A. Pomarol, Phys. Rev. Lett. 73 (1994) 2292; Phys. Rev. D51 (1995) 

6532; D. Matalliotakis and H. P. Nilles, Nucl. Phys. B435 (1995) 115; A. Pomarol and S. 

Dimopoulos, Nucl. Phys. B453 (1995) 83; H. Murayama, Berkeley preprint LBL-36962 

(1995), [hep-ph/9503392| . 
[3] V. Berezinsky, A. Bottino, J. Ellis, N. Fornengo, G. Mignola and S. Scopel, CERN 

preprint CERN-TH 95-206 (to appear in Astroparticle Physics), [hep-ph / 9508249| . 
[4] V. Berezinsky, A. Bottino, J. Ellis, N. Fornengo, G. Mignola and S. Scopel, CERN 

preprint CERN-TH 96-42, hep-ph/ 9603342 . 



[5] J. Ellis, G.L. Fogli and E. Lisi, Phys. Lett. 286B (1992) 85, Nucl. Phys. B393 (1993) 
3, Phys. Lett. 324B (1994) 173 and 333B (1994) 118; P. Langacker and M. Luo, Phys. 
Rev. D44 (1991) 817; M. Carena and C.E.M. Wagner, Nucl. Phys. B452 (1995) 45. 
[6] G. Altarelli, R. Barbieri and F. Caravaglios, Phys. Lett. 314B (1993) 357. 
[7] D. Garcia and J. Sola, Phys. Lett. B354 (1995) 335. 
[8] P. H. Chankowski and S. Pokorski, Phys. Lett. 366B (1996) 188. 
[9] G.L. Kane, R.G. Stuart and J.D. Wells, Phys. Lett. B354 (1995) 350. 
[10] P. H. Chankowski and S. Pokorski, preprint IFT-96/6 ( [hep-ph 96033 10| ). 
[11] M.D. Hildreth, talk at the XXXI Rencontres de Moriond, March 1996. 
[12] A. Djouadi, G. Girardi, C. Verzegnassi, W. Hollik and F.M. Renard, Nucl. Phys. B349 
(1991) 48; M. Boulware and D. Finnell, Phys. Rev. D44 (1991) 2054; A. Blondel and 
C. Verzegnassi, Phys. Lett. B311 (1993) 346. 
[13] J. Erler and P. Langacker, Phys. Rev. D52 (1995) 441. 
[14] L. Roszkowski and M. Shifman, Phys. Rev. D53 (1996) 404. 
[15] M. Shifman, Mod. Phys. Lett. A10 (1995) 605. 
[16] D. Garcia and J. Sola, Phys. Lett. B357 (1995) 349. 

[17] J. Ellis, J.L. Lopez and D.V. Nanopoulos, CERN preprint CERN-TH/95-314, [Eepj 
| ph/951228% 

[18] L. Rolandi, H. Dijkstra, D. Strickland and G. Wilson, representing the ALEPH, DEL- 
PHI, L3 AND OPAL collaborations, Joint Seminar on the First Results from LEP 1.5, 
CERN, December 1995; ALEPH Collaboration, preprint CERN-PPE/96-10; OPAL Col- 
laboration, preprint CERN-PPE/96-019 and CERN-PPE/96-020. 

[19] A. Djouadi, P. Zerwas and J. Zunft, Phys. Lett. B259 (1991) 175. 

[20] J. Wells and G. Kane, Phys. Rev. Lett. 76 (1996) 869. 

[21] We thank Piotr Chankowski for providing us with the outputs of the \ 2 fit of Ref. |§, 

employed in our analysis. 
[22] K. Griest and D. Seckel, Phys. Rev. D 43 (1991) 3191. 

[23] A. Bottino, V. de Alfaro, N. Fornengo, G. Mignola and S. Scopel, Astroparticle Physics 2 
(1994) 77; A. Bottino, V. de Alfaro, N. Fornengo, G. Mignola, S. Scopel, and C. Bacci 
et al. (BRS Collaboration), Phys. Lett. B295 (1992) 330. See also the review paper 
by G. Jungman, M. Kamionkowski and K. Griest, Phys. Rep. 267 (1996) 195, and 
references quoted therein. 



9 



[24] L. Bergstrom and P. Gondolo, preprint OUTP-95-38P, |hep-ph/ 95 10252 



[25] For reviews covering also the experimental aspects, see: J. R. Primack, D. Seckel and 
B. Sadoulet, Annu. Rev. Nucl. Part. Sci. 38 (1988) 751; P. F. Smith and J. D. Lewin, 
Phys. Rep. 187 (1990) 203; R. Bernabei, Riv. N. Cim. 18 (1995) N.5, and L. Mosca, 
invited talk at TAUP 95 (Toledo, September 1995) to appear in Nucl. Phys. B (Proc. 
Suppl.). 

[26] F. J. Kerr, D. Lynden-Bell, Mon. Not. R. Astr. Soc. 221 (1986) 1023. 
[27] P.J.T. Leonard and S. Tremaine, Ap.J. 353 (1990) 486. 

[28] E.I. Gates, G. Gyuk and M. Turner, Fermilab preprint FERMILAB-Pub-95/090-A, 

|hep-ph/9505039| . 
[29] R.A. Flores, Phys. Lett. B215 (1988) 73. 

[30] T.K. Gaisser, G. Steigman and S. Tilav, Phys. Rev. D34 (1986) 2206. 

[31] S. P. Ahlen et al, Phys. Lett. B195 (1987) 603; D. O. Caldwell et al., Phys. Rev. Lett. 

61 (1988) 510. 
[32] D. Reusser et al, Phys. Lett. B255 (1991) 143. 

[33] M. Beck (Heidelberg-Moscow Collaboration), Nucl. Phys. B (Proc. Suppl.) 35 (1994) 

150; M. Beck et al, Phys. Lett. B336 (1994) 141. 
[34] E. Garcia et al, Nucl. Phys. B (Proc. Suppl.) 28 A (1992) 286; M. L. Sarsa et al., Nucl. 

Phys. B (Proc. Suppl.) 35 (1994) 154; E. Garcia et al, Proc. "The Dark Side of the 

Universe", eds. R. Bernabei and C. Tao (World Scientific, Singapore 1994), p. 216. 
[35] A.K. Drukier et al., Nucl. Phys. B(Proc. Suppl.) 28A (1992) 293; I. R. Sagdev, 

A. K. Drukier, D. J. Welsh, A. A. Klimenko, S. B. Osetrov, A. A. Smolnikov, Nucl. 

Phys. B (Proc. Suppl.) 35 (1994) 175. 
[36] We thank Angel Morales for providing us with the listings of the COSME [34] and 

TWIN |55 germanium spectra. 
[37] For a discussion of this point, see Ref. |§] and references quoted therein. 
[38] An analysis, along the general lines of the present paper, based on experimental bounds 

on up-going muon fluxes from the centre of the Earth and from the Sun (indirect neu- 

tralino detection) is in preparation. 



10 



Figure Captions 



Figure 1 - Neutralino relic abundance tt x h 2 versus Rdirect, for tan/? = 1.4. This scatter 
plot displays only supersymmetric configurations that satisfy the condition Qh 2 < 1. The 
rate Rdirect reported in this figure was evaluated, using set I of Table 1. The vertical line 
corresponds to the experimental bound given in Eq. ([!]). 

Figure 2 - Neutralino relic abundance versus the neutralino mass, as a scatter plot over 
all the configurations of the best-fit analysis, for tan/3 = 50. 

Figure 3 - Scatter plots of the rate Rdirect as a function of m x . The rate Rdirect shown 
in section (a) has been evaluated using set I in Table 1, the rate Rdirect shown in section 
(b) using set II. Configurations of different cosmological relevance are denoted by different 
symbols: configurations with a relic abundance Q x h 2 > 0.1 are denoted by circles, those 
with a relic abundance in the range 0.025 < n x h 2 < 0.1 are denoted by crosses, the other 
configurations by points. 

Figure 4 - Plot of Rj, versus x 2 . The solid line delimits the region which contains the 
whole sample of the best-fit configurations of Ref. P|. The dashed curve delimits the region 
where those configurations of the original best-fit sample excluded by R^rect are located. 
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